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We present an analysis of Q,, H, A and (f> spectra from Au+Au collisions at -ysjvjv ~ 200 GeV 
in terms of distributions of effective constituent quarks at hadronization. Consistency in quark ra- 
tios derived from various hadron spectra provides clear evidence for hadron formation dynamics as 
suggested by quark coalescence or recombination models. We argue that the constituent quark distri- 
bution reflects properties of the effective partonic degree of freedom at hadronization. Experimental 
data indicate that strange quarks have a transverse momentum distribution flatter than that of 
up/down quarks consistent with hydrodynamic expansion in partonic phase prior to hadronization. 
Our extracted parton transverse momentum distributions at the hadronization provide a unique 
constraint on the partonic evolution history. After the AMPT model is tuned to reproduce the 
strange and up/down quark distributions, the model can describe the measured spectra of hyperons 
and (j) mesons very well where hadrons are formed through dynamical coalescence. 
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Recent data from the Relativistic Heavy Ion Collider 
(RHIC) at BNL have demonstrated the formation of a 
hot and dense partonic matter in high-energy nuclear 
collisions at RHIC [E S [li 13 • Measurements of nuclear 
modification factor Rcp and elliptic flow V2 for identi- 
fied particles in the intermediate transverse momentum 
(pt) region, 2 < px < 5 GeV/c, exhibit a number of 
constituent quark (NCQ) scaling [1, S S B H- Such 
scaling can be explained by quark coalescence or recom- 
bination models (lol . [ill [l3 | , which provide an intriguing 
framework for hadronization of bulk partonic matter at 
RHIC. The essential degrees of freedom at the hadroniza- 
tion seem to be effective constituent quarks which have 
developed a collective elliptic flow during the partonic 
evolution. The elliptic flow V2 of the constituent quarks 
at hadronization can be characterized by hadron V2 
scaled by the NCQ, W2/NCQ vs. pr/NCQ [1 1 " 
and the hadron elliptic flow results from the sum of con- 
stituent quark collectivity. In this paper, we examine the 
constituent quark number scaling in transverse momen- 
tum spectra and use experimentally measured fi, S, A 
and (p Pt spectra to explore their connections to possible 
quark distributions at hadronization. 

We focus our analysis on fJ, S and (j) spectra in contrast 
to previous coalescence or recombination model calcula- 
tions, e.g. [13, [ll|, where pion and kaon spectra have 
been used to obtain the thermal parton component. The 
final state spectra of common hadrons like pion, kaon 
and protons in nucleus-nucleus collisions represent cumu- 
lative contributions from partonic dynamics, hadroniza- 
tion and kinetic evolution in hadronic stage. Collec- 
tive radial fiow can significantly alter the hadron px dis- 
tribution and ordinary hadrons may decouple from the 
hadronic evolution at a very late stage depending on the 
hadronic interaction cross sections J^,J^]. Multi-strange 
hadrons, <^s, Ss and fls, are predicted to have a relatively 
small hadronic interaction cross section Um, UM ■ These 



hadrons carry the information of partons directly from 
the hadronization stage with little or no distortion due 
to hadronic evolution. In addition, there is no decay feed- 
down contributions to the CI and </> spectra while for ordi- 
nary pion, kaon and protons the majority of the observed 
yield comes from decay production of resonance and weak 
decay states. Multi-strange hadrons offer unique advan- 
tages to probe properties of the partonic degrees of free- 
dom at hadronization. 

Quark coalescence or recombination models have been 
used extensively in explaining RHIC data recently [lOl . 
[ill . [1^ . There are some common features in the in- 
termediate Pt region below 5 GeV/c in these models: 
i) baryons with transverse momentum px are mainly 
formed from quarks with transverse momenta Pt/-^, 
whereas mesons at px are mainly produced from par- 
tons with transverse momenta ~ PtI^- ii) The produc- 
tion probability for a baryon or meson is proportional 
to the product of local parton densities for the con- 
stituent quarks. We argued that the ri(pT/3)/0(pT/2) 
ratio can reflect the strange quark distribution prior to 
hadronization as the 17 baryon consists of three valence 
strange quark (sss) while the meson carries hidden 
strangeness (ss). The S(pt/3)/0(pt/2) ratio will reflect 
the light quark information since the S baryon consists 
of one light valence quark plus two strange quarks. We 
have assumed that the strange and anti-strange quark 
distributions are the same and the particle formation dy- 
namics are dominated by coalescence or recombination 
of bulk partons. Our approach is consistent with the 
recombination model calculation by Hwa and Yang [l^ 
using (5-function approximation for recombination and is 
valid in the intermediate px region up to 5 GeV/c for 
baryons and 4 GeV/c for mesons. Our approach allows 
us to extract properties of partons at hadronization di- 
rectly from recent statistically improved measurements 
of multi-strangeness hadrons at RHIC [1, [^ . The valid- 
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ity of our approach can also be tested by experimentally 
examining parton distributions derived from independent 
ratios of various hadrons. The extracted parton px distri- 
bution provides unique constraints on partonic evolution 
history up to the hadronization epoch. 



FIG. 1: (color online) Left panel: The derived quark pr dis- 
tributions in Au-|~Au collisions at RHIC. Gray bands are fit- 
tings with hydrodynamics inspired model (see text for de- 
tails). Hatched area is the scaled range for the meson 
dN/d(pT/2) divided by the extracted strange quark dN/dpTs 
distributions. Right panel: similar as left panel but plotted 
as KEt scale. 

Fig. □ presents the ratios of ~^M2)^ ^^^^ '^<I>{pt/2) 
as a function of pT/riq, where Uq is the number of con- 
stituent quarks. The data points are from STAR mea- 
surements for Au-f Au collisions at ^snn— 200 GeV and 
the ratios are calculated at the quark- number scaled px- 
The shape of these ratios represent the strange and down 
quark pr distributions. Within the statistical and sys- 
tematic uncertainties there is no significant variation in 
the quark pT distributions at hadronization for collision 
centralities from 0-5% to 40-60%. The collision geome- 
try factor is cancelled out in these ratios. The feed-down 
contribution from higher state resonance S(1530) hasn't 
been included yet and will be discussed in detail in the 
following. 

In order to characterize the quark px distribution, hy- 
drodynamics motivated functions J^] have been used to 
fit the derived quark distributions, permitting extraction 
of model parameters characterizing the bulk freeze-out 
temperature {Tth) and collective radial flow velocities 
{vt)- For practical propose, we follow Fries's model pl| . 
assuming vt to be independent of source radii and az- 
imuth angle with the radial expansion velocity profile 
as assumed in other models [l^. Fittings to the ratios 
derived from hadrons at different centralities simultane- 



ously yielded vt = (0.54 ± 0.13)c and Tth = (131 ± 48) 
MeV for strange quarks with "X^ /ndf — 10.7/12, and vt = 
(0.36 ± 0.19)c and Tth = (170 ± 40) MeV for light quarks 
with x^/ndf = 28.1/41. If we fixed the parameter Tth — 
170 MeV, we obtained vt = (0.43 ± 0.03)c with x^/ndf 
= 11/12 for strange quarks and vt — (0.36 ± 0.02)c with 
/ndf = 26.6/41 for light quarks. Varying the freeze-out 
temperature Tth to 160 (180) MeV , vt values increased 
(decreased) ~ 10% for both strange and light quarks. 
These results are consistent with the picture that strange 
quarks may undergo a stronger hydrodynamical expan- 
sion in partonic phase than light quarks due to the larger 
effective quark mass; e.g., in Fries' model calculation the 
assumed quark masses are 460 MeV for s quarks and 260 
MeV for u/d quarks [11|. A recent proposed transverse 
kinetic energy scaling [17[, KEt = — 'm-O: has been 
applied in this analysis as well (m^ is the hadron trans- 
verse mass while mo is the rest mass). The consistency 
in the KEt scaling, as shown in right panel of Fig. [2 
provides further indication for coalescence or recombina- 
tion of quarks where these quarks must have undergone a 
partonic evolution possibly described by hydrodynamics. 

The physical picture to relate the particles ratios to 
strange and up/down quark distributions can be fur- 
ther tested with independent hadron ratios. In this pic- 
ture ratios of Q{pt/3)/E{pt/3) and E{pT / 3) / A{pt / 3) 
should have a similar shape since both represent the 
ratio of s/d quark distributions. Fig. [5] shows the s/d 
ratios extracted from the fl, S~ and A spectra from 
central Au-|-Au collisions at ^/sWn= 200 GeV, where 
the hadron px has been scaled by the NCQ. The raw 
17(pt/3)/S(pt/3) and E{pt/3)/A{pt/3) ratios have a 
similar pT/rig shape indicating the validity of our ap- 
proach. Recombination model calculation by Fries et 
al. 11] predicted a consistent shape between s/d quark 
ratio and the n{pT/3)/E{pT/3) or E{pt/3)/ A{pt/3) ra- 
tio after removing the different spin degeneracy factor. 
The calculated s/d ratio as a function of quark pT de- 
viates somewhat in shape from our parameterized curve 
based on experimental data. The overall agreement is 
reasonably well because of large experimental uncertain- 
ties involved. There is a normalization offset, presumably 
due to the fact that the model calculation does not in- 
clude all resonance decay contributions. 

Attempt to correct for feed-down contributions re- 
sulted to large statistical errors. The S~ spectrum was 
corrected for feed-down from higher state resonance de- 
cays based on the measured E!°(1530) spectrum [l^. The 
contribution to the S~ spectrum from S'^(1530) decays 
was 46% ± 14%, while the feed-down contribution from 
decays was negligible. The A spectrum has already been 
corrected for S and fl weak decays, while the S contri- 
bution has not been subtracted This is shown as red 
circles in Fig. [21 The E{pT / 3) / A{pT / 3) ratio is about a 
factor of 3 smaller than the 0(pt/3)/2(pt/3) ratio. The 
feed-down contribution from S is significant. The contri- 
bution to the A spectrum from E(1380) decays was 26% 
± 5.9% [l^. The T,^ production has not been measured 
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FIG. 2: (color online) Up panel: The s/d ratio derived from 
n, E~ and A spectra in central Au+Au collisions at ^snn~ 
200 GeV. Red circles represent the s/d ratio without E decays 
correction while gray hatched region is the corresponding ratio 
range after taking into E decays correction. The blue boxes 
represent the uncertainty in H(1530) feed-down estimation. 
Theoretical predictions (see text for details) have also been 
plotted for comparison. Low panel: similar as up panel but 
plotted with KEt scale. 



well experimentally at RHIC. Using thermal model pa- 
rameters fit from central 200 GeV Au-I- Au data [23| , the 
THERMUS thermal model gave a primordial E°/A ratio 
of 0.67 and a final ratio of S° to inclusive A to be 0.36 [2]| . 
On the other hand, string fragmentation model predicted 
a primordial ratio of 1/3 [22]. Assuming the S° feed- 
down is pt independent, its contribution to the A spectra 
was assumed to be 25% - 36%. The E{pt/3)/A{pt/3) 
ratio after taking into account the S decay contribution 
is consistent with the f7(pT/3)/S(pT/3) data [cf. Fig. [2]. 

The large uncertainty due to the feed-down contribu- 
tions will mostly change the normalization scale with- 
out significantly altering the shape of the spectrum in- 
volved in our ratio calculation. The similarity in the ra- 
tios of Q{pt /3) /E{pt /3) and S(pt/3)/A(pt/3) supports 
our picture that one can derive quarks distributions at 
the moment of hadronization from these multi-strange 
hadrons. Another independent check is to use the strange 
quark distribution from f2 to ratios, and derive the anti- 
strange quark distribution from (j) meson data divided by 



the strange quark distribution. The gray hatched area in 
Fig. [1] shows the anti-strange quark distribution, which 
is consistent with the strange quark distribution. 

The consistency in s/d ratios derived from different 
hadron species provides clear evidence for hadron for- 
mation dynamics as suggested by quark coalescence or 
recombination models. The s/d ratios derived from the 
n{pT/3)/E{pT/3) and E{pT / 3) / A{pt / 3) data both in- 
crease with pt for (pt/3) < 1.0 GeV/c, and approach 
saturation for (pt/3) > 1.0 GeV/c. This pt dependence 
may indicate that strange quarks have developed a col- 
lective flow stronger than that of light quarks by the time 
of hadronization |13l . [l3 | . Comparison with Fries's model 
calculation and fit parameters extracted in Fig.[l]are con- 
sistent with stronger radial flow for strange quarks. The 
similarity of the s and d quark KEt distributions [cf. 
bottom panel of Fig. [2])] may indicate a hydrodynamic 
behavior during the partonic evolution as suggested by 
Ref. [l3|. But the exact nature of KEt at the parton 
level is yet to be determined. 

Theoretical models for hadron production in nucleus- 
nucleus collisions at RHIC typically involve initial con- 
ditions arising from parton scatterings, partonic evolu- 
tions, hadronization and hadronic evolutions. Different 
theoretical paradigms separate partonic from hadronic 
processes. Theoretical uncertainties due to hadronization 
scheme and hadronic evolution are major issues for quan- 
titative description of properties of the QCD medium 
created at RHIC. Hydrodynamic calculations, for exam- 
ple [1^, often employs Cooper- Frye [l^l hadronization 
scheme. Fragmentation models or coalescence formation 
models may also be used which lead to different final 
state hadron momentum distribution and azimuthal an- 
gular anisotropy. The hadronic evolution processes have 
been added to hydrodynamic models as an afterburner 
and have been shown to si gnifi cantly alter the spectra 
shapes of ordinary hadrons [2q|. By using effective con- 
stituent quark degrees of freedom and their pT spectra, 
it is hoped that major uncertainties from hadronization 
and hadronic evolutions may be avoided in model calcu- 
lations. We will use a multiphase transport (AMPT) 
model [1^ to illustrate that our derived strange and 
up/down quark distributions, representing a cumulative 
effect from initial conditions through partonic evolution, 
play an important role in determining the final state 
hadron momentum distribution. Details of the AMPT 
model can be found in Ref. [l^ and won't be elaborated 
here. 

The AMPT model has been used extensively to explain 
the RHIC data; e.g., by using parton scattering cross 
sections of 6-10 mb, the AMPT model with string melt- 
ing scenario was able to reproduce both the centrality 
and transverse momentum (below 2 GeV/c) dependence 
of the elliptic flow and pion interferometry measured in 
Au-HAu coUisions at ^51^= 130 GeV at RHIC [13, HI]. 
It can reproduce the measured pt dependence of both 
V2 and V4 of mid-rapidity charged hadrons and V2 for 
0-meson in the same collisions at y^sjvw=200 GeV as 
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FIG. 3: (color online) The AMPT model calculation com- 
pared with the data in central Au+Au collisions at ^sjvjv = 
200 GeV. The data have been scaled by factors for clarity. 



well |29l . |30|. Furthermore, it can produce a conic 
emission pattern through partonic cascading in central 
Au+Au coUisions at RHIC [3l|. These successes prob- 
ably reflect the fact that the AMPT model has the in- 
gredients for general features of partonic space-time evo- 
lution in these collisions. However, the model failed to 
reproduce the hadron transverse momentum spectra at 
RHIC as shown of dashed lines in Fig. [S] This may be 
due to the fact that partons in the AMPT model have 
not undergone the initial partonic evolution stage where 
partons would have developed a large radial flow. We 
remedy this deficiency in the model by modifying the 
original quark distributions before the partonic evolu- 
tion at the moment of HIJING interaction ceased. 
We effectively introduced an initial transverse momen- 
tum distribution to partons from HIJING string melt- 
ing in the AMPT scheme, and let these partons follow 
the partonic evolution as described by Zhang's Parton 
Cascade (ZPC) model We tune the initial parton 

transverse momentum distribution so that after the ZPC 
evolution the parton distributions match our extracted 
strange and up/down quark distributions. The coales- 
cence hadronization scheme in the AMPT model has 
been used for hadron formation. We refer to Ref. [s^l 
for the details of partonic evolution history in the AMPT 
model. 

The introduction of the initial parton transverse mo- 
mentum distribution may be considered as an effective 
way to include early non-equilibrium partonic cascade 



effect, not modelled by the string formation and melting 
scheme. It is also possible that this is an indication of 
insufficient parton cascading cross sections in the ZPC 
model ^] where only pQCD processes have been in- 
cluded. It is important to note that in our empirical 
approach the strange quarks seem to develop a harder 
Pt distribution than up/down quarks during the parton 
evolution. This may require to take into account the ef- 
fective strange quark mass in the evolution process. 

Fig. [3] shows results from the modified AMPT model 
calculation for hyperons and ((> mesons. Our new AMPT 
calculation can reproduce the measured px spectra for 
multi-strange hadrons at mid-rapidity very well. We note 
that the hadronization process in the AMPT model is 
based on the coordinate space information (i.e., two near- 
est quark and antiquark are combined into mesons and 
three quarks or antiquarks are combined into baryons or 
antibaryons that are closest to the invariant masses of 
these parton combinations). In the coalescence scheme 
used by the AMPT model the partons must be close 
in coordinate space and there could be a broad range 
for momentum difference between coalescing partons. 
This is somewhat different from the ones by other 
groups [13, [111, where partons are coalescing in mo- 
mentum space, based on which we have obtained our 
parton transverse momentum distributions empirically. 
The consistency between our new AMPT calculation and 
the data indicates that an essential ingredient is the dis- 
tribution functions of effective constituent quarks which 
readily turn into hadrons in coalescence or recombina- 
tion calculations. Our extracted strange and up/down 
quark distributions provide an unique constraint on the 
partonic evolution history. Theoretical calculations of 
partonic cascading can be compared with our quark dis- 
tributions without involving complicated hadronization 
and hadronic rescattering processes. 

In summary, we have presented constraints on trans- 
verse momentum distributions for the effective con- 
stituent quarks at hadronization of the bulk partonic 
matter produced at RHIC. Our results suggest that 
strange quarks may have developed a collective radial 
flow stronger than that of light quarks during the ini- 
tial partonic evolution. The coalescence model as im- 
plemented in the AMPT model can faithfully reproduce 
the measured multi-strange hadron transverse momen- 
tum spectra at RHIC when our derived quark distribu- 
tions are used at hadronization of the partonic matter. 
The validity of our approach to explore quark transverse 
momentum distributions at hadronization has also been 
tested with independent particle ratios. Our approach in 
complement with the constituent quark number scaling 
in elliptic flow provides a means to measure quantita- 
tive quark properties at hadronization of bulk partonic 
matter. 
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